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ABSTRACT 


The  plasma  sheath  simulation  technique  developed  earlier  has 
been  used  to  investigate  the  effects  of  various  types  of  sheath  disconti¬ 
nuities  and  inhomrgeneities  on  a  slot  antenna  radiation  pattern  and  input 
impedance. 

In  the  case  of  semi-infinite  and  finite-extent  homogeneous 
plasma  sheaths^  the  radiation  pattern  and  the  impedance  characteristic 
are  investigated  with  regard  to  the  geometry  of  the  discontinuity  and 
its  proximity  to  the  slot  antenna.  It  is  foimd  that  the  radiation  pat¬ 
tern  is  imaffected  by  the  discontinuity  as  long  as  the  ratio  of  the  dis¬ 
tance  between  the  source  and  the  discontinuity  to  the  sheath  thickness  is 
of  the  order  of  20  or  greater.  When  the  sheath  is  made  finite  in  extent, 
with  a  discontinuity  on  each  side  of  the  slot,  the  resulting  pattern  is 
seen  to  be  a  superposition  of  the  effects  of  each  discontinuity  acting 
separately.  The  input  impedance  of  the  slot  exhibits  only  small  varia¬ 
tions  when  the  ratio  of  the  discontinuity  separation  to  the  sheath 
thickness  is  of  the  order  of  4  or  greater.  For  smaller  separations, 
the  impedance  exhibits  significant  variations,  strongly  dependent  on  the 
geometry  of  the  discontinuity.  The  effects  produced  by  inhomogeneities 
are  found  to  be  similar  to  those  of  a  homogeneous  sheath  with  a  certain 
equivalent  dielectric  constant. 

For  both  the  homogeneous  and  inhomogeneous  sheaths,  favorable 
comparison  with  the  available  theory  is  obtained. 
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1,  Introduction 


The  recent  interest  in  propagation  of  electromagnetic  waves 
through  a  plasma  sheath  has  generated  a  large  number  of  papers^  mostly 
theoretical  in  nature,  which  solve  for  the  fields  of  various  radiators 
in  the  presence  of  idealized  plasma  sheath  environments.  Most  of  these 
papers  make  assumptions  concerning  the  nature  of  the  discontinuities  and 
inhomogeneities  in  the  plasma  sheath  in  order  to  make  the  problem  more 
tractable  analytically.  For  example,  Tamir  and  Oliner  [j-96^  approxi¬ 
mate  a  discontinuous  plasma  she#‘th  by  terminating  the  sheath  at  a  finite 
distance  from  the  radiator  with  a  termination  perpendicular  to  the  ground 
plane  and  to  the  plasma-air  interface.  As  will  be  seen  later,  this  con¬ 
figuration  allows  an  approximate  analytical  solution  to  be  made  by  means 
of  a  Kirchhoff-Huygens  integration  over  the  discontinuity  plane.  Rusch 
[l963^  and  Tyras  |j.965bJ  formulated  mathematical  models  for  the  dielectric 
constant  in  a  stratified  inhomogenecus  medium,  with  the  latter  obtaining 
a  closed  form  expression  for  the  radiation  fields » 

Little  or  no  experimental  work  in  this  area  has  been  reported, 
primarily  due  to  the  handling  difficulties  and  instabilities  inherent  in 
actu..l  laboratory  plasmas.  Recently,  however,  a  plasma  sheath  environ¬ 
ment  was  successfully  simulated  using  real  dielectric  materials  jlyras 
et  al.,  1965aj  .  The  simulation  method  allowed  the  experimental  determi- 
mation  of  the  radiation  patterns  and  input  admittances  for  the  problem  of 
a  slot  antenna  in  an  infinitely  conducting  plane  clad  with  a  homogeneous, 
lossless,  and  isotropic  plasma  layer  of  infinite  extent. 

This  simulation  method  provides  a  laboratory  geometry  readily 
suitable  for  modification  so  that  discontinuities  and  inhomogeneities  in 
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the  plasma  sheath  may  also  be  simulated.  In  particular,  discontinuities 
of  the  "plate"  tj^je  jlamir  and  Oliner,  196^  as  well  as  inhomogeneities 
composed  of  discrete  vertical  stratifications  can  be  easily  simulated. 

It  is  the  purpose  of  this  report  to  describe  in  detail  these 
simulated  inhomogeneous  and  discontinuous  plasma  sheath  environments  and 
to  present  their  respective  radiation  patterns  and  input  admittances.  In 
the  case  of  the  radiation  patterns,  the  experimental  results  are  compared 
with  the  approximate  analytical  results  of  Tamir  and  Oliner  for  the  dis¬ 
continuous  plasma  sheath  and  with  the  asymptotically  exact  results  ob¬ 
tained  from  a  saddle  point  integration  for  the  inhomogeneous  plasma 
sheath. 


2.  Theoretical  Background 
2,1  Semi-Infinite  Plasraa  Sheath 
A  common  idealized  model  of  the  plasma  sheath  surrounding  a  re¬ 
entry  vehicle  is  seen  in  figure  1,  with  Lj^  and  L2  considered  infinite  in 
length.  It  consists  of  a  perfectly  conducting  plane  covered  with  a  homo¬ 
geneous  plasma  of  infinite  extent,  Tamir  and  Oliner  I19623,  among  others, 
have  analyzed  this  problem  and  have  shown  that  the  field  in  the  plasma 
region  and  the  near  field  in  the  air  region  are  dominated  by  a  single 
leaky  wave  and  that  the  radiation  field  has  a  maximum  at  an  angle  closely 
corresponding  to  the  critical  angle  of  geometric  optics,  as  long  as  the 
plasma  dielectric  constant  remains  positive. 

To  study  the  effect  of  finiteness  on  the  radiation  patterns, 
Tamir  and  Oliner  [l96-^  proposed  the  model  of  figure  1,  with  L2  consid¬ 
ered  infinite,  finite,  and  «  90°,  In  solving  for  the  approximate 
radiation  fields,  it  was  assumed  that  the  fields  incident  on  the  discon” 
tinuity  plane  were  just  those  present  in  the  case  of  the  infinite  plasma 
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slab,  namely,  the  field  due  to  the  leaky  wave.  The  problem  is  then  one 
of  an  aperture  with  a  prescribed  field  distribution.  A  Kirchoff-Huygens 
integration  over  the  fields  on  this  aperture  yields  the  far  field  and 
consequently,  the  radiation  pattern,  R(0.). 


R(9) 


K(e) 
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and  where 


cr  =•  kdcosG  (k  a  free  space  wavenumber) 
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(2) 


i  =  (L/d)cot0  (0-  =  critical  angle) 

kp  s*  k^d  (k^  =  transverse  wavenumber  in  the  air  region) 

kep  **  kze*^  (k^g  =  transverse  wavenumber  in  the  plasma  region) 

and  k  and  k  satisfy  the  transverse  resonance  relation 
p  ep 
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When  the  discontinuity  plane  and  the  source  are  separated  by  a 
large  distance,  the  first  term  in  K(0)  becomes  dominant;  and  the  radiation 
pattern  is  just  that  of  an  infinite  plasma  layer.  Once  the  ratio  of  the  dis¬ 
continuity  separation  to  the  sheath  thickness  is  made  less  than  20,  the  radi¬ 
ation  pattern  begins  to  change  appreciably.  This  change  is  composed  of  a 
broadening  of  the  major  lobe  and  a  shifting  of  the  peak,  all  on  the  disconti¬ 
nuity  side  of  the  source  (0  70).  In  addition,  increased  end-fire  and  broad¬ 
side  radiation  is  observed.  The  exact  nature  of  this  change  is  strongly 
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dependent  upon  the  plasina  sheath  para=:eters  (gp  and  d). 


2.2  InhoEOgeneous  Plasaa  Sheath 
A  possible  sodel  for  the  inbo>20geneity  profile  thought  to  exist 
in  the  plassa  sheath  surrounding  a  reentry  vehicle  is  shovn  in  figure  2a. 
This  profile  can  be  approxlnated,  for  analytical  as  veil  as  experimental 
purposes,  by  constructing  a  sheath  of  discrete  layers  of  homogeneoi*s  die¬ 
lectrics  as  In  figures  2b  and  2c.  As  the  layers  become  finer,  the  dis¬ 
crete  Eiodel  will  approach  the  continuous  model. 

Consider,  then,  an  infinitely  conducting  ground  piane  covered 
with  N  layers  of  homogeneous  dielectric  ssaterial  of  infinite  extent  and 
excited  by  an  infinite  narrow  slot  antenna,  as  in  figure  3.  By  symmetry, 

there  is  no  variation  of  the  fields  with  the  x-coordinate,  or  ~ 

^  2> 

Consequently,  the  only  component  of  the  magnetic  field  is  the  x-ccmpcaeat, 
th 

Then  in  the  a  layer,  satisries  the  source-rree  wave  equation 
2  2 

V  H  +  k  H  =  0  (4) 

n  Xr 


'n 


where  k  =  k  and 

n  o  V  Ti 


a 


-■»  it 

An  e  ~  time  dependence  has 


been  assinaed  and  suppressed  throughout.  Since  the  geometry  is  infinite 
in  the  y-direction,  a  Fourier  transform  pair  on  the  y-coordinate  can  be 
defined  as  follows: 
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e^^da. 
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Thus*  "  =  iQ  and  —  =  -Q  ^  and  equation  (4)  becomes 
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Tbe  solution  of  (7)  is  foisid  to  be 
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Tba  Ajj  sad  3^^  reprsseac  tbe  arolitcdes  of  trzassitted  and  reflected  vsvbSj 
respectively.  In  the  air  regica^  there  is  c-o  reflected  wave;  and  eqszacioa 
(9)  becccses 


(10) 


The  coaosoaeats  of  tbe  electric  field  are  given  by 


2^^ 
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d  y 


(lla) 


(lib) 


The  boundary  conditions  are  the  continuity  of  the  tangential 
E  and  H  fields  at  the  dielectric  interfaces  and  the  prescribed  E  field  at 
the  sotirce.  This  can  be  sunaarized  by 


at  z  =  Q^_i 


(12a) 


Vl  -  ^  -  '^n-l 


(12b) 
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Application  of  (13)  on  equation  (iia)  yields 
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Tfcas  rrv*  cuantiQr  s  (A.  -  E,)  is  recognized  as  ti^  Fourier  transfora  of 

1  i  r 

eic  =  VcCy);  or.  frcca  ecuaticns  (5)  and  (6)  . 
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vbicb  can  be  rewritten  as 
tj=  e  V 


•H  = 


r 

1  1 
1  1  -  (Bi/A,) 


(15  b) 


Using  the  boundary  conditions  for  the  interaediate  interfaces^  equations 
(12a)  and  (i2b)  becoae 

Is  ^  t  “is  ^  “is  _  ftc.  \ 
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Tbs  spplLicatlon  of  eqaations  (i6j  at  tbs  last  Interface,  where  3^*0, 
rescits  in  tbs  foilowiag  expressions: 


Ss-l 


and 
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Equations  (16)  can  now  be  ssanlpulated  into  tbe  following  recursion  re¬ 
lationships  : 

(18) 
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These  recursion  relationships  can  now  be  progra=Ered  on  a  coa- 
puter  to  yield  the  desired  traasaission  coefficient,  Ajj.  The  procedure 
for  calculating  Ajj  uas  the  following: 

1.  and  its  reciprocal  were  calculated  froa  eauatica  (17a) 

Bn-1 

and  were  stored. 

2.  The  result  was  put  into  equation  (IS),  where  a  series  of  itera- 

f  ft 

tions  were  perforaed  to  calculate  Each  CR  and 

^--2  ^  ^ 

its  reciprocal  were  stored. 

3-  Aj^  was  then  calculated  using  equation  <15b;.  This  result  was 
then  placed  in  equation  (19),  where  iterations  were  perforaed 
to  obtain  A2, ....A^_,. 

4.  A-,  was  then  obtained  froa  eauation  (17b). 

Evaluation  of  the  integral  representation  of  the  transaitted 

field  can  now  be  carried  out  by  aaking  the  following  changes  of  the 

variables  in  equation  (7): 

y  =  csin9 

z-d  5=  pcos9 

Ct  =  k  sin^  , 
o 

If  the  last  aediua  is  air,  then  c,.  =  1  and  equation  (iO)  becoaes: 

"lx 


\  = 


/  o 


N  /2f 

le  point  integration  is  used  to  evaluate  this  expression 


The  method  of  saddle  po 
for  the  radiation  field.  The  resulting  expression  is: 


li.  =  e^ -^of"’  &  cos9  A 


■N 


0  =  9 


Since  Ajj  is  ultimately  evaluated  at  <>  =  9,  the  value  of  each  s^ 
the  computer  program  is  just  s^  =  I  €n  -  sin^e  . 


(21) 
used  in 
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3.  Experisentai  Systems  and  Results 
3.1  Sisconcinuous  Plasaa  Sheath 


replaced  the  problem  of  the  plasaa-alr  interface  with 


I- 


The  plasca  sheath  si=uiatioa  systen  described  Tyras  et  al. 

that  of  an 

air-liquid  dielectric  interface.  The  dielectric  used  was  Aroclor  1232,  a 
fluid  with  a  dielectric  constant  of  2.78  and  a  low  loss  tangent.  The  fluid 
was  contained  in  a  seat  -cylindrical  plexiglass  tank  with  a  dielectric  con¬ 
stant  of  2.59.  This  provided  an  essentially  reflection-free  piexiglass- 
aroclor  interface.  It  was  shown  then  that  the  air  layer  between  the 
ground  plana  and  the  tank  was  equivalent  to  a  plasna  sheath  with  a  die¬ 
lectric  constant  of  0.36,  and  that  the  Aroclor  region  was  equivalent  to 
the  air  region  above  the  plasna  sheath.  Thus  the  radiation  patterns  ob¬ 
tained  for  the  experimental  system  were  expected  to  exhibit  the  same 
functional  relationships  as  those  obtained  fxoa  the  analytical  solution 
of  the  actual  problem. 

That  this  was  a  feasible  method  for  a  plasma  sheath  simulation 
was  shown  by  the  excellent  agreement  between  the  theoretical  and  experi¬ 
mental  radiation  patterns.  The  variation  of  the  input  admittance  with 
the  sheath  thickness  was  also  investigated,  and  good  agreement  with  the 
available  theoretical  data  was  again  observed. 

The  extension  of  this  method  to  simulate  finite  and  semi-infinite 
plasma  sheath  environments  is  straightforward,  A  plexiglass  sheet  was 
placed  in  the  space  between  the  tank  and  the  ground  plane  in  order  to 
simulate  the  rectangular  semi-infinite  plasma  sheath  of  figure  1,  with 
considered  Infinite  and  equal  to  90°,  Radiation  patterns  were  then 
obtained  for  various  values  of  the  discontinuity  separation  L^.  Similar 
measurements  were  also  made  for  angular  discontinuities  with  equal  to 
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45°  and  15^.  It  is  expected  that  this  type  of  discontlnuivy  is  a  laore 
realistic  isodel  of  the  environments  encountered  in  reentry. 

Figures  4a  through  4f  snow  the  experimental  and  theoretical 
radiation  patterns  obtained  for  the  rectangular  discontinuity.  The  fea¬ 
tures  of  major  importance  in  these  patterns  are  the  following: 

1.  The  radiation  pattern  is  relatively  unaffected  as  long  as 
the  ratio  of  the  discontinuity  separation  to  the  sheath 
thickness  is  at  least  20* 

2.  As  the  discontinuity  plane  is  moved  closer  to  the  source^ 
the  major  lobe  of  the  radiation  pattern  on  the  disconti¬ 
nuity  side  of  the  source  (0^0)  is  seen  to  broaden  and  the 
critical  angle  is  seen  to  shift.  In  addition,  increased 
broadside  and  end-fire  radiation  is  observed. 

3.  While  it  is  not  shown  on  all  the  figures,  the  radiation 
pattern  on  the  other  side  of  the  source  (9^0)  is  affected 
very  little  by  the  presence  of  the  discontinuity. 

These  features  are  completely  in  line  with  the  predictions  of  Tamir  and 
Oliner  Il964]  on  Che  basis  of  their  leal^-wave  analysis. 

In  comparing  the  theoretical  and  experimental  radiation  patterns, 
it  is  to  be  noted  that  no  outstanding  point-by-point  agreement  was  reached. 
This  is  due  in  part  to  the  inherent  errors  arising  from  the  approximate 
leal^-wave  method  of  analysis  for  obtaining  the  radiation  fields.  This 
error  is  usually  apparent  at  broadside  for  certain  combinations  of  plasma 
dielectric  constant,  Cp,  and  sheath  thickness,  d,  which  is  readily  apparent 
in  figure  4b,  However,  the  general  qualitative  nature  of  the  effects  of 
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the  discontinuity  on  the  radiation  patterns  is  in  good  agreenent  with  the 
theory  in  view  of  the  three  features  noted  above. 

With  a  view  toward  siisulating  rsore  realistic  dis continuities , 
the  geoaetry  of  figure  1,  with  finite  L  and  L,,  was  also  investigated. 

It  was  observed  that  the  effects  of  the  angular  discontinuities  were  of 
the  sane  general  type  as  those  of  the  rectangular  discontinuity.  Again^ 
no  appreciable  changes  were  noticed  for  separations  greater  than  10  wave¬ 
lengths;  and  the  changes  noted  for  soaller  separations  were  of  the  sane 
type  as  those  observed  previously.  Figures  5a  and  5b  show  the  radiation 
patterns  obtained  for  an  angular  discontinuity  with  ii^finite  and  = 

45°.  The  effect  of  the  addition  of  a  second  discontinuity  to  sinulate 
the  environnent  of  figure  1  is  seen  in  figures  6a  id  6b.  It  is  evident 
that  the  resulting  radiation  pattern  is  essentially  a  superposition  of 
the  effects  of  each  discontinuity  acting  separately. 

The  effects  of  discontinuities  on  the  Icput  admittance  were  also 
investigated.  Measurements  were  made  at  the  input  to  the  tapered  waveguide 
section  of  the  antenna,  exactly  as  reported  earlier 

The  resulting  admittance  characteristics  are  seen  in  figures  7a  and  7b. 
Figure  7a  shows  the  variation  of  input  admittance  with  the  discontinuity 
separation  with  only  one  discontinuity  present.  It  is  seen  that  the  con¬ 
ductance  remains  essentially  constant  if  the  ratio  of  the  discontinuity 
separation  to  the  sheath  thickness  is  at  least  4,  while  the  susceptance 
remains  constant  over  the  whole  range  of  the  variation  in  Lj^.  The  effect 
of  the  second  discontinuity  on  the  input  admittance  is  seen  in  figure  7b. 
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3.2  Inhonogeneous  Plasma  Sheath 
As  mentioned  in  Section  2.2^  the  discrete-layered  model  for  the 
inhomogeneous  plasma  sheath  is  a  convenient  one  analytically  as  well  as 
experimentally.  The  profiles  shown  in  figures  2b  and  2c  are  easily  simu¬ 
lated  in  the  laboratory  using  the  earlier  scaling  technique  jlyras  et  al., 
1965^ -  Table  I  below  shows  the  properties  of  the  various  materials  used 
to  simulate  these  profiles. 


Table  1 

(Data  for  25°C  and  3x10^-1x10^  cps) 

Material  Actual  e  Scaled  g 


Styrofoam 

1,00 

0.360 

Teflon 

2.08 

0.748 

Polyetlq^lene 

2.25 

0.813 

Polystyrene 

2.54 

0.913 

Plexiglass 

2.59 

0.935 

Figures  8a  and  8b  show  the  theoretical  and  experimental  radiation 
patterns  obtained  for  these  two  profiles.  In  figure  8a,  good  agreement 
with  theory  is  observed,  with  the  main  discrepancy  occurring  in  the 
slightly  different  critical  angles.  In  figure  8b,  a  difference  of  2  db 
is  noted  in  the  attenuation  at  broadside,  while  there  is  excellent  agree¬ 
ment  at  the  critical  angle.  In  each  of  these  figures,  the  most  signifi¬ 
cant  point  to  observe  is  the  general  shape  of  the  radiation  patterns.  It 
seems  to  indicate  that  an  inhomogeneous  plasma  sheath  of  the  type  shown  in 
figures  2b  and  2c  acts  effectively  like  a  homogeneous  plasma  sheath  with 
some  average  dielectric  constant.  This  average  value  of  gp  can  be  readily 
computed  from  the  relation,  sin  0^  =  gp,  where  0^  is  the  critical  angle  ob¬ 
tained  from  the  radiation  pattern.  This  should  hold  at  least  as  long  as 
the  minimum  value  of  the  dielectric  constant  in  the  sheath  remains  positive. 
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The  effects  of  discontinuities  in  the  inhomogeneous  plasma  sheath 
were  also  investigated.  The  configurations  of  figures  9a  and  9b  were  simu¬ 
lated  in  the  laboratory  using  the  materials  of  Table  I.  The  radiation  pat¬ 
terns  for  various  values  of  the  discontinuity  separation,  Lj^,  are  seen  in 
figure  10,  These  patterns  exhibit  the  same  characteristics  as  those  of 
the  discontinuous  homogeneous  sheath. 

The  input  admittance  of  the  antenna  in  the  presence  of  these 
semi-infinite,  inhomogeneous  environments  was  also  obtained.  Figure  11 
shows  the  conductance  and  susceptance  as  a  function  of  the  discontinuity 
separation.  As  noted  before  in  the  case  of  the  homogeneous  sheath,  the 
variation  is  very  slight  for  separation-to-thickness  ratios  of  4  and 
larger.  It  is  only  when  the  discontinuity  is  situated  closer  to  the 
source  that  significant  variations  are  observed. 


4.  Conclusions 

The  extension  of  the  earlier  results  j^ras  et  al.,  1965a^  to 
include  more  realistic  radiating  system  environments  has  succeeded  in 
predicting  the  qualitative  effects  of  the  discontinuities  and  inhomogene¬ 
ities  in  the  plasma  sheath  on  the  radiation  patterns  and  input  admittances. 
In  the  case  of  the  finite  and  semi-infinite  homogeneous  plasma  sheath,  it 
has  been  shown  that  the  analytically  approximate  method  of  Tamir  and  Oliner 
[1964]  is  a  valid  one  to  use  in  order  to  predict  the  essential  features  of 
the  radiation  pattern.  It  was  observed  that  the  radiation  p  tttern  was 
relatively  unaffected  if  the  discontinuity  separation-to-sheath  thickness 
ratio  was  about  20  or  greater.  For  smaller  separations,  the  radiation  pat¬ 
tern  was  observed  to  broaden  and  the  major  lobe  to  shift.  In  the  case  of 
a  stratified  inhomogeneous  sheath,  this  report  has  shown  that  the 
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Fig.  1.  Plasma  Clad  Slot  Antenna  Geometry 
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Fig.  1,  Plasma  Clad  Slot  Antenna  Geometry 
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Fig.  2  Models  of  Plasma  Sheath  Inhomogeneity  Profile:  (a)  Contiraous 

Model^  (b;  Four  Layers  Approximation^  (c)  Six  Layers  Approximation 
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Fig.  3.  Infinite  Inhomogeneous  Plasma  Sheath  Geometry 
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Fig.  4d  Experimental  and  Theoretical  Radiation  Patterns  for  Semi-Infinite 

Homogeneous  Plasma  Sheath 
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Fig.  4e  Experimental  and  Theoretical  Radiation  Patterns  for  Seml-Inflnlt 

Homogeneous  Plasma  Sheath 
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Fig.  5a  Experimental  Radiation  Pattern  for  Semi-Infinite  Homogeneous 

Plasma  Sheath 
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Fig.  5b  Experimental  Radiation  Pattern  for  Serai-Infinite  Homogeneous 

Plasma  Sheath 
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Fig.  8b  Experimental  and  Theoretical  Radiation  Patterns  for  Infinite 
Inhomogeneous  Plasma  Sheath  (6  Layers) 
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Fig,  9a  Semi-Infinite  Inhomogeneous  Plasma  Sheath  Geometries 
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Fig.  9b  Semi-Infinite  Inhomogeneous  Plasma  Sheath  Geometries 
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Fig.  11  Admittance  Characteristic  for  Semi-Infinite  Inhomogeneous  Plasma  Sheath 
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